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Abstract 

We discuss effects induced by beyond leading order contributions to photon pair production. 
We point out that next to leading order contributions to the fragmentation component of the 
signal lead to a change in the shape of distributions. This is already mildly visible in present 
Tevatron data though stringent isolation criteria tend to suppress the fragmentation component 
considerably. We expect the effect to be experimentally confirmed in future data samples with 
higher statistics which would serve as a precision test for QCD. 
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1 Introduction 



The production of prompt photon pairs with large invariant mass is the object of continuing exper- 
imental studies ranging from fixed target £jL 0, El, to colliders energies [EL M, 0], especially by the 
CDF and DO experiments at the Tevatron. The word "prompt" means that these photons do not 
come from the decays of hadrons such as ir° , i], etc. at large transverse momentum. Prompt photons 
may be produced according to two possible mechanisms: either they take part directly to the hard 
subprocess (direct mechanism), or they result from the fragmentation of large transverse momen- 
tum partons (fragmentation mechanism). For a general discussion of the production mechanisms of 
prompt photon pairs, see Ref. jjj. Besides offering an interesting probe of the short distance QCD 
dynamics, di-photon hadroproduction is the irreducible background for the search of a neutral Higgs 
boson in the channel H — > 77 in the mass range 90 - 140 GeV at the LHC, which provides a strong 
motivation for its extensive study. 

Collider experiments at the Tevatron and the forthcoming LHC do not perform inclusive photon 
measurements. Indeed, the inclusive production rates of jir and 7r°7r° pairs with large invariant 
mass are orders of magnitude larger than the one of prompt photon pairs, to which they provide the 
background. Thus the experimental selection of prompt photons requires the use of isolation cuts. 
The isolation criterion used by the Tevatron experiments CDF Q and DO || is schematically the 
following. A photon is said to be isolated if, in a cone in rapidity and azimuthal angle about the 
photon direction, the amount of deposited hadronic transverse energy E T ad is smaller than some 
value Et max fixed by the experiment: 

E T ad <E Tmax inside {y - y 7 ) 2 + (0 - < R 2 (1) 

The actual experimental criteria, which are defined at the detector level, are much more complicated 
than the simple one given by Eqn. (|l|), see |5|, jjj. They cannot be implemented in a partonic 
calculation. The schematical criterion of Eqn. (0) is a modelization of their effects at the parton 
level. The topic of isolation of photons based on the above criterion is extensively discussed in 
the literature ||, [| |l^, [ll| [l^, [l3|. Besides the rejection of the background of secondary photons, 
the isolation cut also reduces the prompt photons from fragmentation. The stringent isolation 
requirements used by CDF and DO suppress the fragmentation component quite a lot; therefore, 
the effects of the fragmentation component are generally neglected beyond their lowest order in 
perturbation theory. Yet they can affect various observables in a sizeable way. For example, in 
the case of single photon production, it yields a significant increase in the transverse momentum 
distribution in the lower range of the px spectrum |ll], [l^, . 

In a recent work, we presented a full next to leading order (NLO) study of di-photon hadropro- 
duction based on a computer code of partonic event generator type, DIPHOX [Q. In this short 
article, relying on this general study, we illustrate the impact of beyond leading order contributions, 
especially to the fragmentation component, in di-photon production at the Tevatron. 



2 Impact of higher order corrections on di-photon spectra 

Higher order QCD contributions to a given observable lead to a refinement of theoretical predic- 
tions in several aspects. If already the Born level contains the strong coupling a s , higher orders 
stabilize dependencies on arbitrary scales stemming from the truncation of the perturbative series. 
In the hadronic production of direct photon pairs this is not the case though it is in general true 
for the fragmentation parts. Furthermore, the inclusion of the NLO contributions allow for new 
initial states. The additional production channels are typically the reason for large K factors. But 
higher order contributions also modify final state observables, as extra particles in the final state 
generate kinematical configurations which are forbidden at leading order and lead to sizeable ef- 
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fectsQ, especially if new collinear situations in the final states are possible. Such effects can be seen 
in observables sensitive to kinematical configurations where the photons are emitted in rather close 
directions to each other. Indeed, such configurations contribute only beyond leading order (at lowest 
order, transverse momentum conservation forces the photons to be back to back in the transverse 
plane). Moreover, the contribution of such configurations coming from the fragmentation mechanism 
is enhanced, as will be discussed below. 

A typical representative of this familly of observables is the distribution da/dqx of transverse 
momentum of di-photons (qr = ||pr(7i) + Pt(j2)\\)- Fig- @ shows the comparision between the DO 
preliminary data and our theoretical computation taking into account the kinematic and isolation 
cuts of the DO Tevatron experiment. DO requires |y(7i,2)| < 1.0 for the photon rapidities, and the 
lower experimental pr cuts are 

Pt(h) > 14 GeV, p T ( 72 ) > 13 GeV (2) 

which roughly correspond to effective cuts at respectively 

PTmUll) = 14-90 GeV, PTrmn(l2) = 13.85 GeV (3) 

in a partonic calculation^. Furthermore, a lower cut in the two-photon acollinearity is required: 

(2/(71) -y(7 2 )) 2 + ^ 7 >R 2 min (4) 

where </> 77 is the photon-photon azimuthal angle and R m in — 0.3. The isolation criterion is modelized 
according to Eqn. (|l|) with Exmax — 2 GeV, R = 0.4. 

One can see a shoulder at about qr ~ 35 GeV in the distribution measured experimentally. At 
first glance, one may question its significance, given the limited statistics of the Run I data and 
the systematic uncertainties which affect these preliminary data. However a corresponding shoulder 
appears in the theoretical calculation performed with the same binning as the data too. A carefull 
study reveals that this shoulder is not a mere statistical fluctuation of the partonic Monte-Carlo 
generator, nor a binning artefact, but it is a physical effect. 

In order to see this, we split the two-photon phase space into regions, according to the photon- 
photon azimuthal angle 77 : 

7T 7T 

Region I: < </> 77 < — ; Region II: — < 4> 11 < ir. (5) 

Since at leading order, the two body kinematics forces the photons to be back to back in the 
transverse momentum plane (</> 77 = 7r), the leading order contributions to da/dqx come entirely 
from region II, and vanish in region I. With respect to the above artificial slicing (^), a "new channel" 
opens at NLO - and beyond - which populates region I. However, as qr — ||pt(7i) + Pt(72)||, the 
asymmetric cuts on the pr's of the photons imply that the contribution to da/dqx from region I 
vanishes for q T < q Tm in with q Tm in = (Pt mm (7i) + Pr mm (72)) 1/2 , so that this "new channel" 
opens only above qTmin — 20.34 GeV. It is instructive to examine separately the direct, and the 
fragmentation components^] above this threshold. 

2 Higher order corrections can also reveal new infrared sensitive issues, in particular inside the physical spectrum, 
as it is the case for the qx distribution at qx = Ex max |t| . 

3 The experimental cuts at measured values 14 and 13 GeV respectivelly in the DO data are not corrected for the 
electromagnetic calorimeter absolute energy scale. See H for more details. 

4 We call direct the mechanism by which both photons take part to the hard subprocess, one fragmentation the 
one by which one photon takes part to the hard subprocess whereas the other comes from fragmentation, and two 
fragmentation the one by which both photons come from fragmentation. We call fragmentation the sum of these last 
two components. Actually, the severe isolation cut used by DO suppresses the two fragmentation component down to 
a completely negligible level. Consequently, we will ignore it in the present discussion. 
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The direct component is displayed on Fig. ||(a). As qr increases above qrmin, the size of the 
opening slice of available phase space belonging to region I becomes comparable with the one of 
region II. The direct component, which starts at NLO for qx > qrmim is by definition free of any 
final state collinear singularity. It thus involves only smooth variations with (2/(71) — 2/(72)) and 
</> 77 . Therefore, after a transition regime extending up to q^ < 35 — 40 GeV, the magnitudes of 
the contributions from regions I and II become comparable. Moreover the R m in dependence, which 
comes from region I, is weak, 0(R^ nin ) at small R m in- In summary the approximate doubling of 
available phase space due to the opening of region I above qx min roughly amounts to a doubling of 
the direct component. 

On the contrary, the fragmentation component is collinearly enhanced when the two photons 
come close to each other. Indeed, without the cut ([|), contributions such as the one pictured in 
Fig. |^ would make the fragmentation component logarithmically divergent for all qx larger than 
qr Um = PTmin(li) + Pr minilz) with q T u m ~ 28.75 GeV. This divergence would be produced when 
the photon considered as "direct" would become collinear to the one from fragmentation^. The 
acollinearity cut (|J) forbids this configuration; yet a collinear enhancement survives in region I as 
a memory of this regime. As the slice of available phase space in region I opens, the fragmentation 
contribution from this region rises initially faster than the contribution from region II decreases 
above qx min ■ This generates a peak in the fragmentation component about qr ~ 30-35 GeV, untill 
the contribution from region I decreases too. As a reflect of the collinear enhancement, the smaller 
the value of R m in , the higher the peak. For qr min above 25-30 GeV, the fragmentation term comes 
mostly from region I, contrarily to the direct case. This is illustrated on Fig. ^| which compares the 
case Rmin = 0.3 experimentally used and what would happen for R m in = 0.01, an artifical value 
chosen on purpose to emphasize this effect. 

Due to the fact that the fragmentation component is severely suppressed by the isolation cut (|l]) 
with respect to the direct one, and that the value R m in = 0.3 experimentally used leads to a rather 
moderate enhancement of this fragmentation component, the effect of the latter is limited to a 10% 
modulation effect in the total (= direct + fragmentation) distribution, so that the shoulder is mainly 
due to the direct component. This can be seen in Fig. [| showing the ratio total/direct. However, 
the effect would be more spectacular with lower values of R m in- For illustrative purpose, Fig. [| also 
shows what would be the situation with the value R rn in = 0.01. 

The collinear enhancement in the fragmentation component affects also the distribution da / d<j) 11 
of azimuthal angle between the two photons in a pair, in the lower range of the <j) 11 spectrum, cf. 
Fig. ^. The contribution of the direct component decreases monotonically while remaining finite 
when </> 77 < 7r decreases towards 0. Instead, the fragmentation contribution to this observable would 
diverge logarithmically at </> 77 = without the acollinearity cut (Q), as can be inferred from Fig. ^. 
The fragmentation contribution "feels" this collinear divergence by rising again when </> 77 decreases 
below ~ 1 towards R m in- This divergence is forbidden by the constraint (jij). The latter produces a 
turnover when 77 < R m im which can be seen on the right plot on the right of Fig. |6| in the lower 
range of the spectrum. 

One may wonder about the significance of these effects compared with the theoretical scale 
uncertainties associated with the arbitrariness of the renormalization and factorization scales used 

5 From a theoretical point of view, the subtraction of the corresponding double collinear divergence for the {partonic 
emitter + two photons} system in the partonic calculation would require a new type of fragmentation functions of 
partons into "twin collinear photons", extending the procedure corresponding to the case of fragmentation into a 
single photon. Such collinear photons would actually be detected as single photon events. Most probaby they would 
contribute to the single photon production rate at a negligeable level, at least in the Tevatron energy range, due to 
a relative suppression factor oc 0(a em ) with respect to the typical single photon production rate. In practice, the 
experimental cuts on the acollinearity (3/(71) — J/(72)) 2 + </> 77 of the two photons, or on the invariant mass of the pair, 
prevent from having to consider this regime. 
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in the partonic calculation^. Concerning the qx distribution, the direct contribution to the tail 
of the distribution above qx min starts at NLO (the Born contribution to the direct component is 
concentrated at qx = by transverse momentum conservation) . As for the c/> 77 distribution, both the 
direct and fragmentation contributions to the tail of the distribution away from </> 77 = n start at NLO 
too (the respective Born terms being concentrated at 77 = it). For each of these two observables, 
NLO is the effective lowest order. Consequently, it is fair to say that these uncertainties are rather 
large. However scale uncertainties act roughly as an overall effect: they do not generate local 
distortions comparable to the pattern of effects studied above. We insist that the latter are physical 
effects, which could in principle be unambiguously observed among data with higher statistics and 
an improved understanding of systematics, as it may hopefully be the case in the forthcoming Run 
II of the Tevatron. 

3 Conclusions and perspectives. 

In this note we have emphasized NLO effects, coming in particular from the fragmentation contri- 
bution to di-photon production. We have pointed out that a new collinear sensitive situation arises 
in the fragmentation part at next to leading order, namely when the two photons become more 
and more close to each other. This leads to an enhancement of the respective signal. Stringent 
isolation criteria reduce the fragmentation contribution considerably and acollinearity cuts between 
the two photons forbid collinear photons in the experimental situation. Consequently these next to 
leading order effects from fragmentation are not very significant in the present data. However we 
have shown that the account of these contributions is required for an accurate description of the 
shapes of di-photon observables. They thus may become relevant for a precise understanding of the 
high statistic data to be collected during the forthcoming Run II at the Tevatron. 
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Figure 1: qr distribution of photon pairs. Black dots: DO data white diamonds: average values 
of the NLO calculation (DIPHOX code) in the corresponding experimental bins. The curve is a 
spline interpolation between the theoretical average. 
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Figure 2: Origin of the qr shoulder in the theoretical calculation. Plot (a) shows the direct component 
(open squares) split into the phase space regions 77 < tt/2 (full triangles) and ci 77 > it/ 2 (open 
triangles) for the experimentally used value R m in — 0.3. Plot (b) is the same as (a) but for R m in = 
0.01 to show the sensitivity of the effect to this collinear cut. Plots (c) and (d) show the corresponding 
histograms for the fragmentation component where the enhancement due to photon collinearity is 
clearly visible. 
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Figure 3: Kinematical configuration for which the fragmentation contribution is collinearly enhanced 
when the two photons become close to each other. 
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Figure 4: Top: the qr spectrum (open squares) split into direct (full triangles) and fragmentation 
part (inversed triangles) for R m in — 0.3 (a) and R m in = 0.01 (b). Bottom: the ratio of the total qr 
distribution divided by the direct one for Rmin = 0.3 (c) and R m i n = 0.01 (d). 
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Figure 5: 77 distribution of photon pairs. Black dots: DO data J^/, white diamonds: average values 
of the NLO calculation (DIPHOX code) in the corresponding experimental bins. The curve is a 
spline interpolation between the theoretical average. 
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Figure 6: Top: i/ie theoretical 77 distribution (open squares) split into direct (full triangles) and 
fragmentation part (inversed triangles) for R m in = 0.3 (a) and R m in = 0.01 f&j. Bottom: the ratio 
of the total </> 77 distribution divided by the direct one for R m in = 0.3 fc,) and Rmin = 0.01 f'rfj. The 
collinear enhancement of the fragmentation part for small </> 77 is clearly visible. 
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